Structural and hyperfine characterization of cr-phase Fe-Mo alloys 
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A series of nine samples of cr-Feioo-^Moa; with 44 < x < 57 were synthesized by a sintering 
method. The samples were investigated experimentally and theoretically. Using X-ray diffraction 
techniques structural parameters such as lattice constants, atomic positions within the unit cell and 
populations of atoms over five different sublattices were determined. An information on charge- 
densities and electric field gradients at particular lattice sites was obtained by application of the 
Korringa-Kohn-Rostoker (KKR) method for electronic structure calculations. Hyperfine quantities 
calculated with KKR were successfully applied to analyze Mossbauer spectra measured at room 
temperature. 

PACS numbers: 33.45.+X, 61.05.cp, 61.43.-j, 71.20.-b, 71.20.Be, 71.23.-k, 74.20.Pq, 75.50.Bb, 76.80.+y 



I. INTRODUCTION 

A cr-phase (space group P4 2 /ninm, 30 atoms per unit 
cell), which belongs to a class of the so-called Frank- 
Kasper phases 1 - characterized by a high coordination 
numbers (CN=12-15, for a) can be formed in alloy sys- 
tems by a solid state reaction at elevated temperatures. 
About 50 examples of a were up-to-now found in such 
binary systems, including the Fe-Mo one. A possibility 
of occurrence of the cr-phase in Fe-Mo was first revealed 
in 1949 by Goldschmidt 2 . Studying a ternary system of 
Fe-Cr-Mo he had identified a new phase stable in the 
temperature interval of 1180 - 1540 °C. The phase had 
some similarities to the already-known cr-phase in Fe-Cr. 
Its crystallographic structure and parameters were def- 
initely identified five years later—. Wilson was the first 
who determined the lattice parameters and distribution 
of atoms over the five sublattices for cr-FesoMoso sample 4 . 
Further investigations depicting the structure and sites 
occupancy by Fe and Mo atoms were continued by other 
investigators, too^. Heijwegen determined borders for 
the existence of the cr-phase in the Fe-Mo system^. Fi- 
naly, the phase diagram which is regarded as the most 
actual was proposed by Guillermert-. 

In the Fe-Mo system, the cr-phase is one of three Frank- 
Kasper phases that can occur in that system. The other 
two ones are represented by \i (16 atoms per unit cell), 
and R (54 atoms per unit cell). They have different com- 
positions and can be formed in different ranges of temper- 
ature than cr— In the available literature, there is little 
information relevant to physical properties of the cr-phase 
in Fe-Mo. Even the data on its crystallographic struc- 
ture in the whole concentration range of its occurrence is 
not complete. In these circumstances systematic studies 
aimed at filling this gap are justified, all the more so such 
studies have been recently carried out for this phase in 
Fe-Cr and Fe-V systems 1 ^. A comparison of correspond- 
ing data obtained for different binary alloy systems con- 
taining iron, Fe-X, where the cr-phase can be formed, is 
of interest per se, and, additionally, it may contribute 
to a better understanding of its formation mechanism. 



The konowledge of the letter is of a great importance in 
the light of a detrimental effect of a cr-phase presence in 
stainless steels used in various branches of industry as 
valuable construction materials. In this paper results of 
a systematic study of the lattice parameters, sites occu- 
pancies, Fe-site charge-densities and electric field gradi- 
ents obtained using various experimental tools viz. X-ray 
diffraction (XRD) and Mossbauer spectroscopy (MS), as 
well as electronic structure calculations (KKR method) 
are presented in the whole concentration range of the 
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FIG. 1: (Online color) Parts of of selected fitted X-ray diffrac- 
tograms recorded at 294K on the cr-phase sample of (a) 
Fe55Mo45 and (b) Fe45Mos5. The solid line stays for the best- 
fit obtained with the procedure described in the text. Peak 
positions for n and a (a) as well as for pure cr-phase (b) are 
indicated. A difference diffractograms are shown, too. 
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FIG. 2: (Online color) Dependence of the lattice parameters 
a (circles, blue) and c (diamonds, red) on Mo-content, x, as 
determined from the X-ray diffractograms recorded at 294K. 



cr-phase existence in the investigated alloys. 



II. EXPERIMENTAL 

It was not possible to obtain samples in a pure er-FeMo 
alloys by using the procedure that was previously suc- 
cessfully applied in the case of Fe-Cr and Fe-V alloys i.e. 
an isothermal annealing of ingots obtained by melting 
constituents 11 . The Fe-Mo ingots treated in that way 
were always two-phase, the second phase being /j. Suc- 
cessful was, however, a procedure described by Bergman 3 
viz. a sintering. For that purpose, powders of elemental 
Fe (99.9at%) and Mo (99.95at%) were mixed together in 
adequate proportions. Next, 2 g tablets for each com- 
position were fabricated by pressing the mixtures. The 
tablets were subsequently isothermally annealed at 1430° 
C during 6 hours, and afterwords quenched into liquid 
nitrogen. For XRD and Mossbauer spectroscopic mea- 
surements, both carried out at room temperature (RT = 
294K), the pellets were attrited into powder in an agate 
mortar. Altogether nine samples with the following nom- 
inal concentration of Mo: 44, 45, 47, 48, 49, 51, 53, 55 
and 57 at% were prepared. Based on the XRD patterns, 
an example of which is shown in Fig. 1, it was found that 
all samples except two with the lowest Mo content, and 
one with the highest content of Mo were single-phase, 
namely a. The samples containing nominally 44 and 45 
at% Mo had some admixture of /i-phase (Fig. la), and 
the most Mo-concentrated sample contained ~ 4.2wt% 
of undissolved Mo-rich, bcc Fe-Mo phase. In the latter 
case the content of Mo in the cr-phase was estimated as 
equal to 56.7 at%. 



III. RESULTS 

The powder XRD patterns were collected at RT with 
a D5000 Siemens diffractometer (using Cu K-a radiation 
and a graphite secondary monochromator) from 10° to 
140° in steps of 0.02° in 28. Data were analyzed by 
the Rietveld method as implemented in the FULLPROF 
program^ with 22 free parameters: 6 of them were re- 
lated to a background and positions of lines, 10 param- 
eters represented sites occupancies and atomic positions 
in the unit cell. Remaining 6 parameters described line 
widths, lattice constants and Debye- Waller factors. The 
analysis yielded values of the lattice constants, a and c 
(Fig. 2, Table 1), atomic positions (Table 2), as well as 
occupancies of the sublattices (Fig. 3). Concerning the 
lattice constants, as clearly evidenced in Fig. 2, they 
exhibit a linear dependence on the Mo content, x. Yet, 
the ratio c/a — 0.52233(2) is rather :r- independent. For 
comparison, Bergman found the value of 0.5237 for this 
ratioi The values of the lattice constants given in the 
literatur e) 3 ' 4 ' 13 are in line with ours. This can be taken 
as indication that the nominal compositions of our sam- 
ples can be regarded as very close to the real ones. Also 
the atomic positions found in our analysis agree within 
error limits with those known from the literature 3 -*^. As 
far as the occupancy of the lattice sites is concerned, the 
analysis of the XRD patterns gave evidence that sites A 
and D were exclusively populated by Fe atoms, which is 
slightly different than in the case of the cr-phase in Fe- 
Cr and Fe-V systems where the population of Fe atoms 
at these sites was lying within 80-90% ranged. The re- 
maining three sites were found to be occupied by both 
Fe and Mo atoms, the former being in minority. How- 
ever, the actual occupancy is concentration dependent 
and it decreases linearly with x at a rate characteristic of 
a given site. These results agree qualitatively with those 
reported earliep 4 ^. 




FIG. 3: (Online color) Probability of finding Fe atoms at dif- 
ferent lattice sites in the er-Feioo-iMox compounds, P, versus 
Mo concentration, x. Solid lines stay for the linear fits to the 
data 
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FIG. 4: (a) 57 Fe Mossbauer spectra recorded on a series of 
(T-Feioo-rrMoa; samples at 294K and labeled with the corre- 
sponding i-values. The solid lines are the best-fit to the ex- 
perimental data. The derived isomer shift distribution curves, 
shown in the same sequence as the spectra, are indicated in 
(b). 



Fe-site Mossbauer spectra, shown in Fig. 4a, were 
recorded in a transmission geometry using a standard 
equipment and Co/Rh source for the 14.4 keV gamma- 
radiation. No traces of magnetism at RT were found in 
all the spectra recorded at RT. Despite it is known from 
the XRD experiments that Fe atoms are present at all five 
sublattices, which means that each spectrum should be 
composed of at least five subspectra. However, an unique 
analysis of the measured spectra in terms of five compo- 
nents that could be associated with the five sublattices 
is not possible due to a lack of a well-defined structure 
of them. Instead, one can analyze them in terms of an 
isomer shift distribution, ISD, as explained elsewhere^. 
The ISD-curves obtained as a result of such analysis are 
presented in Fig. 4b. By integration of these curves, 
average values of the isomer shift, (IS), were obtained 
(Table 1). The dependence of the latter on the Mo con- 
tent, x, is illustrated in Fig. 5. It is clear than an increase 
of x results in a linear decrease of (IS) for all samples, 
except for those with the lowest a;-values. The departure 
in the latter case is likely due to the fact that these two 
samples had some small admixture of the phase. The 
decrease of (IS) (increase of the Fe-site charge-density) 
with x is rather weak, as it is equal to 0.0017 mm/s per 
at%. For comparison, in the case of the cr-FeV alloys, the 
rate of (IS) decrease was almost two times higher™. 

A lack of a well-resolved structure of the spectra did 
not allow to decompose them uniquely into the subspec- 
tra corresponding to the five sublattices, hence to deter- 
mine charge-densities and electric field gradients charac- 



TABLE I: Lattice constants (in A), a and c as well as the 
average isomer shift, (IS), (in mm/s, relative to the Co/Rh 
source), for all investigated er-Feioo-ajMoa; samples. 



X 


a 


c 


(IS) 


45. 


9.1598(2) 


4.7844(1) 


-0.242 


46. 


9.1719(2) 


4.7905(1) 


-0.235 


47. 


9.1789(1) 


4.7953(1) 


-0.234 


48. 


9.1866(2) 


4.7987(1) 


-0.236 


49. 


9.1913(1) 


4.8016(1) 


-0.236 


51. 


9.2084(1) 


4.8101(1) 


-0.240 


53. 


9.2251(1) 


4.8188(1) 


-0.242 


55. 


9.2402(1) 


4.8275(1) 


-0.248 


56.7 


9.2489(1) 


4.8322(1) 


-0.252 



teristic of these sublattices. However, as demonstrated 
elsewher o 16 ' 17 , a knowledge of these two spectral quanti- 
ties can be obtained with the help of electronic structure 
calculations. In the present case, the KKR calculations 
were carried out for 23 unit cells with different configu- 
rations of Fe and Mo atoms on the sublattices A, B, C, 
D and E, keeping experimental lattice constants in all 
cases. The configurations were chosen according to the 
following two criteria: (1) the probability of finding an 
Fe atom on a given site was as close as possible to the one 
found experimentally for er-Fe56Mo44, and (2) each pos- 
sible number of Fe atoms being the nearest-neighbors for 
a given lattice site, NNp e , has been taken into account. 

The charge and spin self-consistent Korringa-Kohn- 
Rostoker Greens function method^r— was here used to 
calculate the electronic structure of the Fe-Mo er-phase. 
The crystal potential of the muffin-tin form was con- 
structed within the LDA framework, using the Barth- 
Hedin formula for the exchange-correlation part. The 
experimental values of lattice constants and atomic po- 
sitions were used here. For fully converged crystal po- 
tentials electronic density of states (DOS), total, site- 
composed and /-decomposed DOS (with l max = 2 for 
both types of atoms) were derived. Fully converged re- 
sults were obtained for ~ 120 special k-points grid in the 
irreducible part of Brillouin zone. DOS were computed 
using the tetrahedron k space integration technique and 
~ 700 small tetrahedrons^. More details can be found 



TABLE II: Atomic crystallographic positions, x, y, z, for the 
five lattice sites of the Fe-Mo cr-phase. Since the values were 
weakly composition-dependent, the average values over all 
samples are given. 



Site 


Wyckoff symbol 


X 


y 


z 


A 


2i 











B 


4f 


0.3982(5) 


0.3982(5) 





C 


8i 


0.4635(4) 


0.1297(5) 





D 


8i 


0.7456(9) 


0.0686(15) 





E 


8j 


0.1819(3) 


0.1819(3) 


0.2469(11) 
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FIG. 5: (Online color) The average isomer shift (relative to 
the Co/Rh source), (IS), versus Mo concentration, x, as mea- 
sured (circles) and as calculated (lines). The dashed line 
corresponds to the results obtained with approximation (1), 
while the solid line with approach (2). 



elsewher e) 16 ' 17 , 

The charge-densities, p e , obtained as a result of these 
calculations were analyzed as a function of the number of 
the nearest- neighbor Fe atoms, NNp e . Graphical illus- 
trations of the p e = f(NNpe) dependences are presented 
in Fig. 6. Two types of analysis were carried out. First it 
was assumed the p e -NNp e dependence, was linear, which 
is best seen for the sites D and E. Taking into account 
probabilities of each atomic configuration considered, the 
(JS')-values were calculated for each sublattice, and, fi- 
nally, for each spectrum. The dependence of (IS) on x 
obtained within this approach is shown in Fig. 5 as a 
dashed line. Obviously, it does not agree with the values 
derived from the I S D-curves represented by circles in 
Fig. 5 which means the the linear dependence assumed 
in this approach was not correct. In the second approach, 
the average isomer shift for each subspectrum in a given 
spectrum was calculated as an arithmetic average over 
all atomic configurations considered in the calculations. 
The average IS for the whole spectrum was computed as 
a weighted average over all five sublattices (the weight of 
each subspectrum being its relative spectral area). The 
(IS)(x) dependence determined in this way is shown in 
Fig. 5 as a full line. It is clear that it agrees very well 
with the experimentally found values. 

Based on the measured sublattice occupancies, cal- 
culated potentials and using the extended point charge 
model as outlined elsehwerei 7 -, it was also possible to de- 
termine quadrupole splittings values, QS. As shown in 
Fig. 7, the QS-values are characteristic of a given lat- 
tice site, and for all of them but A, they exhibit a weak 
linear increase with x (for A there is a weak linear de- 
crease). The smallest QS 1 - value of 0.21 mm/s was found 
for the sublattice D, while the largest one, 0.47 mm/s, 
for the site E. The values of QS for the other three sites 
have intermediate values. Noteworthy, similar values and 



relations were found for the cr-phase in Fe-Cr and Fe-V 
system s 16 ' 17 which is reasonable in the light of the same 
crystallographis structure. 

The physical quantities determined in this study viz. 
the site occupancies, the charge-densities (isomer shifts) 
and the quadrupole splittings are sufficient to carry out 
a proper analysis of the Mossbauer spectra in terms of 
the five subspectra corresponding to the five sublattices. 
Only five free parameters were needed to successfully an- 
alyze the spectra. Four of them i.e. a background, spec- 
tral area, line width (common for all five subspectra) and 
the isomer shift of one of the subspectra depend on the 
conditions of the spectra measurements, so they cannot 
be calculated. The fifth free parameter was a propor- 
tionality constant between the calculated electric field 
gradient component, V zz , and the corresponding spectral 
parameter viz. QS. Examples of the spectra modeled in 
this way together with the calculated subspectra corre- 
sponding to the particular lattice sites, are illustrated in 
Fig. 8. A very good quality of the fits achieved can be 
taken as a proof that using the electronic structure calcu- 
lations (KKR-method) the determined spectral parame- 
ters adequately represent values of the hyperfine param- 
eters characteristic of the investigated samples. In other 
words, a combination of the Mossbauer-effect measure- 
ments with the calculations of the electronic structure 
permited to propertly analyze the former in terms of the 
spectral parameters corresponding to the different lattice 
sites. 



IV. CONCLUSIONS 

The following conclusions can be drawn based on the 
results reported in this paper: 

1. The pure cr-phase in the Fe-Mo alloy system can be 
obtained between about 46 and 56.7 at% Mo. 

2. Its lattice parameters a and c linearly increase with 
the Mo content, x. 

3. The c/a ratio is concentration independent and 
equal to 0.52233(2). 

4. The sites A and D are 100% occupied by Fe atoms, 
while the sites B, C and E are populated by both types 
of atoms. The population of Fe atoms on the latter three 
sites is in minority and it linearly decreases with x. 

5. The average isomer shift weakly decreases with x. 
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FIG. 6: (Online color) Fe-site charge density, p e , for five crys- 
tallographic sites versus NNfc- Dashed lines correspond to 
the best linear fit to the data, whereas solid lines represent 
average values of p e 
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6. The calculated charge-densities are characteristic of 
a given sublatticc; the highest one being at site D, and 
the lowest one at site B. 

7. The calculated quadrupole splittings are character- 
istic of a given sublattice: the largest one being at site 
E, and the smallest one at site D. 
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FIG. 8: (Online color) 57 Fe-site Mossbauer spectra recorded 
at 294K on two studied samples viz. with (a) x — 48 and 
(b) x = 55. The best-fit spectrum and five subspectra are 
indicated by solid lines. 



